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Summary 

'TThe  present  study  deals  with  optimal  design  of  arrays 
of  protective  structures.  The  structures,  intended  for 
storage  of  explosive  materials,  consist  of  rectangular 
reinforced  concrete  (RC)  plates,  beams,  and  doors. 

The  design  variables  include:  a)  The  elements1  cross 
sectional  dimensions;  b)  The  structural  configuration, 
c)  The  glometric  location  of  the  structure.  The  cons¬ 
trains  are  related  to  safety  distances,  functional 
requirements,  and  structural  behavior.  The  objective 
function  represents  the  cost,  including  cost  of  mater¬ 
ials,  real  estate,  subgrade,  maintenance,  etc. 

The  main  difficulticj  involved  in  this  problem  stem 
from  the  complex  analysis  and  the  nature  of  the  design 
variables,  constraints,  and  objective  function.  The 
various  types  of  variables  are  of  fundamentally  differ¬ 
ent  nature,  and  some  cf  the  variables  are  discrete. 

The  objective  function  is  neither  differentiable  nor 
continuous. 

App. iximate  behavior  models  are  employed  in  order  to 
simplify  the  analysis.  Since  it  is  not  practical  to 
optimise  all  the  design  variables  simultaneously,  it 
is  proposed  to  use  a  multilevel  optimisation  procedure. 
The  variables  are  optimized  in  different  levels, 
according  to  their  type  anu  nature.  Considerations  for 
choosing  the  levels  are  discussed  and  numerical  examples 
illustrate  the  approach  and  its  practicality. 


1.  Introduction 


Due  to  functional  requirements  the  configurational 
variables  are  discrete  and  the  objective  function, 
representing  the  cost,  is  neither  differentiable  nor 
continuous . 

To  simplify  the  analysis,  approximate  behavior  models 
are  employed.  These  include  empirical  loading  express¬ 
ions,  idealized  constitutive  laws,  and  simplified 
dynamic  models.  A  multilevel  optimization  procedure, 
in  which  the  different  types  of  variables  are  treated 
separately,  is  proposed.  Several  authors,  e.g.  [6], 
proposed  separate  design  spaces  for  geometry  and  cross 
sections  in  truss  optimal  design.  In  general,  such  a 
solution  may  be  viewed  as  a  multilevel  approach  [2,7,8] 
In  the  formulation  presented  in  this  study,  the  cross 
sectional  dimensions  are  optimized  in  the  first  level 
for  a  given  configuration  and  location  of  the  struct¬ 
ure.  As  a  result,  the  structural  elements  can  be 
optimized  independently  in  a  simple  manner.  Data 
banks  of  optimal  elements  are  introduced,  to  be  used 
in  the  higher  levels  of  the  optimization.  J-n  the 
second  level  the  structural  configuration  is  optimized 
for  a  given  location  of  the  structure.  For  each 
candidate  geometry  the  optimal  cross  sectional  dimen¬ 
sions  are  chosen  from  the  data  banks.  The  geometric 
location  of  the  structure  is  selected  in  the  third 
level.  For  each  selected  location,  both  the  struc¬ 
tural  configuration  and  the  cross-sections  are  opti- 
ini  zed . 

The  proposed  approach  combines  efficient  suboptimiza¬ 
tion  for  cross  sectional  variables,  reduction  in  the 
number  of  design,  variables  optimized  simultaneously, 
and  improved  convergence. 


Protective  structure:  intended  to  resist  the  effects 
of  explosions,are  usually  expensive  due  to  the  nature 
of  the  loadings  and  design  criteria.  Therefore, 
improved  designs  may  lead  to  considerable  savings  in 
the  total  cost  of  the  structure.  While  much  work  has 
been  done  in  the  last  two  decades  on  optimum  structural 
design  [1,2],  most  applications  are  limited  to  cross 
sections  optimization.  Examples  of  previous  work  on 
optimal  design  of  protective  structures  include:  cross 
sections  design  [31  or  optimization  of  geometry  and 
cross  sections  [4]  of  RC  slabs  subjected  to  impulse 
loadings;  and  optical  [average]  strength  of  magazine 
doors  [5], 

The  main  difficulties  involved  in  design  of  protective 
structures  stem  from  the  complex  analysis  and  the  nature 
of  the  design  variables,  constraints,  and  objective 
function.  In  general,  the  nonlinear  dynamic  analysis 
must  be  repeated  many  times  and  the  various  types  of 
variables  involved  are  often  of  fundamentally  differ¬ 
ent  nature,  from  both  the  mathematical  and  the  physical 
points  of  view. 

The  present  study  deals  with  optimal  design  of  arrays 
of  protective  structures,  intended  for  storage  of 
explosive  materials.  The  structure  consists  of  RC 
rectangular  plates,  and  the  design  variables  include: 

[a]  The  elements  cross  sectional  dimensions  [plate 
thicknesses  and  amounts  of  steel]; 

[b]  The  structural  configuration  [length  and  width]; 

[c]  The  geometric  location  of  the  structure 
(distances  between  adjacent  structures]. 


2.  ? rob 1 em  Statement 


Consider  the  array  of  RC  rectangular  magazines  shown 
in  Fig,  1.  The  object  is  to  minimize  the  cost  of 
storing  a  unit  weight  of  explosive  material.  It  is 
assumed  that  the  number  of  structures  is  large  and 
all  the  structures  are  identical;  i.e.,  the  optimal 
design  represents  a  standard  magazine.  The  pre  • 
assigned  parameters  include:  materials  properties., 
unit  costs,  height  of  the  structure,  and  spacings 
between  beams.  Uniform  cross  sectional  dimensions  of 
the  RC  elements  and  standard  steel  doors  have  been 
assumed.  The  optimal  design  problem  can  be  stated  as 
follows:  find  the  vector  of  design  variables 
{X}  *  {D},{L},[TC},{AS}  where 


{D}  =  {.D ^ }  (distance  between  adjacent  structures) 


(L)  =  {l  ,L)  (length  and  width  of  the  structure, 

1  1  CZ  K  I  CZ  A  \ 


Lj  £  €  A,) 


{TC}  =  (TCJ(...TCj}  (thickness  of  concrete  plates) 
{AS}  =  {AS  , ...ASj}  (amount  of  steel  in  plates) 


(1) 


such  that 

J 

C  *  \  C./W  -*■  min. 

j  =  l  J 

{DL}  <  {0}  <  {DU} 


(objective  function) 
(bounds  on  distances) 


(2) 

(3) 
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{LL}  <  U)  <  {LU} 
{TCL}  <  {TC}  <  {TCU} 
{pL}  <  {p}  <  tpU> 

{V}  <  {VU} 

{0}  <  [euf 

{oL}  <  to}  <  (aUt 

txL}  <  tx>  <  {tU} 


(bounds  on  structural 
geometry) 

(bounds  on  concrete 
thickness) 

(bounds  on  steel 

percentage) 

(deflection  constraints) 

(rotation  constraints) 

(bending  stress 

constraints) 

(shear  stress 

constraints) 


(4) 

(5) 

(6) 

(7) 

(8) 
(9) 

(10) 


In  the  above  formulation 

Aj,A^  =  sets  of  allowable  discrete  values  of  L.  and  L„, 
respectively  i 

I  =  number  of  elements  in  the  structure 
L,U  =  superscripts  denoting  lower  and  upper  bounds, 
respectively 

C  «  total  cost  of  a  magazine  per  unit  stored 
explosive  material 

W  =  quantity  of  explosive  material  in  a  magazine 
C-  =  the  j-th  component  of  the  objective  function 
J  =  number  of  components  of  the  objective  function 
p  =  steel  percentage 
V  =  deflection 
6  -  rotation  at  the  supports 

a  -  bending  stress 
T  -  shear  stress 

It  should  be  notea  that  (DL)  is  a  function  of  W.  The 
behavior  constraints  (Eqs.  (7)  to  (10))  arc  implicit 
functions  of  the  design  variables,  given  by  the 
analysis  equations. 

Two  types  of  structures,  earth  covered  and  uncovered 
(Fig.  2)  have  been  considered.  The  cost  function 
(Eq.  (2))  includes:  cost  of  materials  (concrete,  steel, 
doors),  real  estate,  subgrade,  maintenance,  etc.  Due 
to  functional  requirements  the  variables  { L}  are  dis¬ 
crete.  The  behavior  constraints  restrict  the  amount  of 
damage  due  to  possible  explsions  in  adjacent  structures. 


3.  The  Analysis  Model 

In  the  discussion  that  follows  only  blast  loadings  have 
been  considered.  It  is  assumed  that  other  effects 
(such  as  fragments)  are  secondary  and  may  oe  checked 
after  the  optimization  process.  Modifications  in  the 
design  can  then  be  made,  if  necessary.  An  approximate 
analysis  model  has  been  employed  with  the  following 
features. 

a)  Loadings 

The  explosive  materials  are  represented  by  an  equi¬ 
valent  TMT  charge  at  the  center  of  gravity  of  the 
charge  [9,10].  The  blast  loadings  due  to  a  possible 
explosion  are  computed  by  the  methods  described  in 
TM  5-1300  [11].  Ths  model  is  based  on  experimental 
results  and  idealized  (such  as  piecewise- linear) 
pressure-time  relations. 

b)  Constitutive  equations  and  resistance  deflection 
functions 

Idealized  elaso-plastic  behavior  models  have  been 
employed.  Also,  it  is  assumed  that  cracking  (but  not 
spalling)  may  occur  in  the  concrete.  The  piecewise 
linear  resistance-deflection  functions  fer  the  various 


elements  are  detemined  by  the  yield  line  theory,  con¬ 
sidering  dynamic  values  of  material  constants  [11], 

c)  Dynamic  response  of  elements 

To  compute  displacements  and  stresses  in  the  elements 
the  medium  is  discretized,  resulting  in  systems  of 
lumpou  masses  and  nonlinear  springs.  A  single  degree 
of  freedom  system  is  assumed  for  each  RC  element,  with 
the  following  equation  of  motion  [11] 


F-p  ■  =  mea  (11) 

in  which  F  *  time  dependent  external  force;  F  *  inter¬ 
nal  force  (resistance);  m  =  mass  of  the  element; 
a  =  acceleration;  =  a  factor  relating  the  value 

of  the  actual  mass  m  and  the  equivalent  mass  me.  The 
value  of  X  is  determined  from  the  principal  mode 
of  vibration?  For  example,  consider  the  sector  of  a 
plate  shown  in  Fig.  3.  The  equation  of  motion 
(rotation  about  the  support)  is 

I 

Fc  -  (EM^  ♦  EMp)  =  -f  a  (12) 


where  c  -  distance  of  the  resultant  force  from  the 
support;  M^,  Mp  =  negative  and  positive  internal 

moments,  respectively;  1  =  moment  of  inertic ;  i  = 

width  of  the  element.  From  Eqs.  (11),  (12)  we  obtain 
I 


K 


m 

Lm  cfcm 


(15) 


For  an  element  consisting  of  a  number  of  sectors 
E(I  /<*) 

K  -  in 

Lm  Em 


(14) 


For  earth  covered  elements  the  corresponding  mass  of 
the  cover  is  also  considered  in  the  equation  of  motion 
[12].  Mult j -degree-of-freedom  systems  may  be  consider¬ 
ed  for  each  element  if  better  approximations  are  re¬ 
quired.  Numerical  algorithms  (such  as  Runge-Kutta)  or 
available  results  111]  are  used  to  solve  the  equations 
of  motion  and  to  evaluate  displacements  and  stresses. 


d)  Computational  considerations 

Analysis  of  a  single  element,  including  calculation 
of  the  loadings  and  solution  of  the  equations  of  motion, 
involves  much  computational  effort  (up  to  ID  sec  Cl  on 
IBM  370/168  were  reported  for  coraple..  loading  histories). 
There  are  several  elements  in  a  structure  and  the 
analysis  usually  must  be  repeated  many  times  during 
optimization.  Therefore,  it  is  proposed  to  reduce  the 
number  of  analyses  in  the  solution  process  by  introd¬ 
ucing  data  banks  of  preoptimized  elements  for  sequences 
of  given  loadings  and  element  configurations. 


4 .  Multilevel  Optimization 

Design  of  a  large  complex  system  usually  involves 
decomposition  into  a  number  of  smaller  subsystems, 
each  with  its  own  goals  and  constraints  [2].  In 
general,  an  integrated  problem  cannot  be  decomposed 
into  subproblems  which  can  be  independently  optimized. 
There  may  be  many  different  ways  of  transforming  an 
optimization  problem  into  a  multilevel  problem.  In 
the  model  coordination  approach  [2,7,8]  used  in  this 
study,  we  choose  certain  variables  called  coordinating 
or  interaction  variables  to  control  the  lower  level 
systems.  For  fixed  values  of  the  coordinating  vari¬ 
ables  the  lower  level  problems  often  become  independ¬ 
ent  and  simple  to  optimize.  The  task  in  the  higher 
levels  is  to  choose  the  coordinating  variables  in  such 
a  way  that  the  independent  lower  level  solutions  are 
optimal.  The  term  "model  coordination"  derives  from 
the  circumstance  that  a  constraint  is  added  to  the 
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problem  in  the  form  of  certain  fixed  interaction  vari¬ 
ables. 

Define  the  vector  of  cross  sectional  design  variables 
(Q)  as 

(Q)  =  (TC), {AS}  CIS) 

and  the  vector  of  configurational  and  geometric  vari¬ 
ables  {R}  by 

{R}  =  {D},{L}  (16} 

It  has  been  noted  earlier  in  section  2  that  the  vector 

of  design  variables  { X}  can  be  partitioned  (Eq.  (1)) 

{X}  =  {R},{Q}  (17) 

In  this  formulation  {R}  is  the  subvector  of  coordina¬ 
ting  variables  between  the  subsystems  and  {Q}  is  the 
vector  of  subsystem  variables,  in  turn  partitioned  as 
follows 

{Q}  =  {Q1),...,{3i) . {Qj}  (18) 

The  subvector  {Q^}  represents  the  cross  sectional 
variables  associated  with  Jie  i-th  element  (subsystem) 
and  1  is  the  number  of  elements.  With  these  defini¬ 
tions,  the  objective  function  of  Eq.  (2)  and  the 
constraints  of  Eqs.  (3)  to  (10)  can  be  expressed  in 
the  general  form 


I 


C  »  F({X})  =  1  F  ({R},{Q  }) 

i=l  1  1 

(19) 

«.  ({X})  =  g.  ({R})  <  0 

A.  f  — 

(20) 

hj({X})  =  h.({R},{Q,};  <  0 

(21) 

in  which  gg  are  constraints  on  the  {R}  variables 
(i.e.,  Eqs.  (3), (4))  and  hp  are  constraints  assoc¬ 
iated  with  the  i-th  element  (Eqs.  (5)  to  (10)).  That 
is,  the  variables  { K>  may  appear  in  all  expressions, 
while  the  variables  CQ^ V  appear  only  in  the  constraints 
associated  with  the  i-th  element,  and  in  the  correspond¬ 
ing  term  of  the  objective  function.  Specifically,  it 
is  assumed  that  the  cross  sectional  variables  of  a 
given  element  affect  only  the  constraints  and  object¬ 
ive  function  component  of  that  element . 

Tile  general  optimization  problem  can  now  be  formulated 
ss  the  following  two  level  problem. 

First-level  problem.  Determine  a  fixed  value 
through  the  constraints 

for  {R} 

{R}  -•  {R°} 

(22) 

Then  the  integrated  problem  can  be  decomposed  into  the 
following  I  independent  first-level  subproblems: 
find  {Qi}  (i=l,...,I)  such  that 

Ci  =  Fp{lRcJ  .{QpJ)  min- 

(23) 

hiC{R°},{Qi})  s  0 

(24) 

Second- level  problem.  The  task  in  the  second 
problem  is  to  find  0'.°}  such  that 

I 

H({R°})  =  l  H. ({R°})  +  min. 
i=] 

level 

(25) 

gk«k°})  1  0 

(26) 

where  IT({R0})  is  defined  by 

kk({R°})  ■=  min  Ch 

(27) 

An  additional  constraint  on  {R°}  is  that  the  first- 
level  problem  has  a  solution,  i.e.,  that  H({ro}) 
exists . 

The  two-level  problem  is  solved  iteratively  as  follows: 

1.  Choose  an  initial  value  for  the  coordinating 
variables  {R°}. 

2.  For  a  given  {R0}  solve  the  1  independent  first- 
level  problems. 

3.  Modify  the  value  of  { R°>  so  that  H({R0})  is 
reduced. 

4.  Repeat  steps  2  and  3  until  min  H({R°})  is  achieved. 

If  all  intermediate  values  for  {r},{q}  are  feasible, 
the  iteration  can  be  terminated  always  with  a  feasible- 
even  though  nonoptimal-soiution,  whatever  the  number  of 
cycles.  This  is  advantageous  from  an  engineering  point 
of  view  and  may  considerably  reduce  the  computational 
effort,  particularly  if  the  object  is  to  achieve  a 
practical  optimum  rather  than  the  theoretical  one. 

Since  the  second  level  variables  ({D}  and  {L},  see 
Eq.  (16))  are  of  fundamentally  different  nature,  it  is 
proposed  to  decompose  the  second-level  problem  into 
two-levels,  such  that  only  the  configurational  varia¬ 
bles  { L)  are  optimized  in  the  second  level,  while  the 
geometric  location  variables  {D}  are  treated  in  a  new 
third-level  problem.  The  three  problems  are  solved 
iteratively  until  the  optimum  is  achieved.  Note  that 
the  I  first-level  subproblems  remain  unchanged 
(Eqs.  (22)  to  (24)).  The  modified  second  and  third- 
levei  problems  are  formulated  as  follows  (Fig.  4). 
Second-level  problem.  For  a  given  geometric  location 


{D)  =  {D°} 

(28) 

find  {Lcj  such  that 

C({D0},{L°},{Q})  -*•  min 

(29) 

gk({D°},{L°})  <  0 

(30) 

Third-level  problem.  Find  { D° }  such  that 

C({D°},(L},{Q>)  •+■  min 

(31) 

gk({D°},{L))  <  0 

(32) 

The  proposed  solution  procedure  is  possible  since  the 
system  by  its  very  nature  c  in  be  decomposed.  The 
loadings  depend  only  on  the  f R}  variables,  therefore 
the  first  level  problems  can  be  solved  for  fixed 
loadings.  The  main  advantage  is  that  I  independent 
simple  subproblems,  are  obtained  in  this  level.  It  has 
been  noted  that  { D‘J>  is  a  function  of  W  ,  which  in 
turn  is  a  function  of  { L} .  Thus,  Eqs.  (3)  are  the 
only  explicit  constraints  which  are  functions  of  both 
•CD}  and  {L}.  The  main  reason  for  choosing  {l}  as  the 
second-level  variables  is  their  discrete  nature.  For 
any  given  {D°)  only  a  limited  number  of  f L)  values 
must  be  considered. The  solution  process  is  shown  in 
Fig.  5. 

5 .  Optimization  Methods  for  the  Various  Levels 

First-level.  For  the  relevant  ranges  of  loadings  and 
element  dimensions,  data  banks  of  preoptimized  elements 
have  been  prepared.  As  a  first  step  a  sequence  of 
loadings,  with  fixed  peak  pressure  and  va-'iablc  dura¬ 
tion,  is  computed  for  a  given  element  conl'i  uration. 

The  cross  sectional  dimensions  are  optimi-eu  by  a 
uirect  search  technique.  For  each  selected  concrete 
thickness  the  amount  of  reinforcing  steel  is  minimized, 
and  the  optimal  thickness  is  evaluated  by  quadratic 
interpolation.  The  data  banks,  to  be  used  in  the 
higher  level  optimization,  are  obtained  by  repeating 
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this  process  for  several  element  configurations. 
Intermediate  solutions  nay  be  evaluated  by  inter¬ 
polation. 

Second-level.  The  optimisation  method  in  this  level 
is  based  on  a  direct  search  in  the  space  of  the  dis¬ 
crete  variables  { L} .  For  each  assumed  Li  value  the 
optimal  L2  is  selected.  Lj  is  then  modified,  with 
the  current  optimal  Lj  (or  the  nearest  feasible  value) 
chosen  as  the  initial  design.  The  Lj  values  are  modi¬ 
fied  until  the  optimum  is  reached. 

Third-level .  Powell's  direct  search  method  [ 13]  is 
applied  in  this  level  with  the  convergence  criterion 

|C({X  })  -  CUX  ^I/Cltx^})  <  ec  (33) 

where  ec  is  a  predetermined  parameter  and  q  is  the 
iteration  number  in  the  third  level.  An  additional 
criterion 

ilV  '  lV}|  <e* 

ensures  that  the  iteration  is  terminated  only  if  the 
condition  of  Eq.  (34)  holds  for  two  successive  itera¬ 
tions. 


6.  Numerical  Examples 

Two  types  of  magazine  have  been  considered  [14] : 
earth  covered  structure  with  a  standard  cover  thickness 
and  an  uncovered  structure.  Fig.  6  shows  some  reasi- 
hle  combinations  of  the  (I.)  variables  (a  =  standard 
spacing,  b  =  required  spacing  to  satisfy  functional 
requirements,  c  =  required  spacing  for  doors).  Fig.  7 
shows  a  typical  design  space  for  the  cross  sectional 
variables,  in  which  As  and  A£  arc  the  amounts  of  re¬ 
inforcing  steel  required  for  tensile  and  compressive 
forces,  respectively.  A  direct  search  in  the  space 
of  d  (the  effective  depth  of  the  cross  section),  as 
shown  in  Fig.  b,  provides  the  optimal  element  for  the 
given  loading  and  configuration.  This  procedure  is 
repeated  for  all  relevant  loadings  and  element  confi¬ 
gurations.  Typical  optimal  costs  for  a  sequence  of 
triangular  pulses  are  shown  in  Fig.  9,  in  which  p  is 
the  peak  pressure  and  the  time  axis  denotes  the 
duration. 

Since  the  second  level  variables  are  discrete,  the 
objective  function  in  the  third  level  is  not  continuous 
(Fig.  10 J .  It  can  be  noted  that  local  optimum  points 
may  exist  and  the  search  is  terminated  only  after 
checking  the  objective  function  in  a  number  of  points, 
to  ensure  that  the  true  optimum  is  achieved.  For  each 
point  in  the  (D)  variables  space  an  optimal  structure 
is  selected.  For  larger  values  of  (D)  the  feasible 
region  in  the  (I.)  variables  space  is  increased,  and 
the  result  might  be  a  change  in  the  discrete  optimal 
(I.)  values.  The  computational  effort  in  the  second 
level  optimization  may  be  reduced  considerably  by 
interactive  design.  Convergence  display  programs  may 
be  used  for  this  purpose.  Fig.  11  demonstrates  the 
relative  cost  against  the  number  of  calculated  (R) 
vectors  (Eq.  (16)).  The  latter  number  is  equal  to  the 
numher  of  accesses  to  the  data  banks.  SeaTch  in(cVsp<ie« 
is  sKovo-n  in  Fig.  12.  Different  init ial  desi gns  or 
directions  in  thft  space  provided  an  identical  opti¬ 
mum.  Earth  covered  structures  were  optimized  in  a 
similar  manner,  with  convergence  demonstrated  in  Fig, 

13.  This  type  of  structure  was  found  to  be  cheaper 
than  the  uncovered  structure. 

The  possibility  of  choosing  the  structural  configura¬ 
tion  (y  as  the  third-level  variables  was  also  checked. 
While  an  identical  optimum  was  obtained,  the  conver¬ 


gence  rate  was  much  slower;  the  number  of  iterations 
was  about  50%  larger.  This  can  be  explained  by  t he 
difficulties  encountered  in  searching  the  optimum  in 
the  discrete  variables  space.  Choosing  11.}  as  the 
second-level  variables,  the  number  of  cheeked  points 
for  a  given  (U)  value  is  limited  hy  the  constraints 
of  Eq.  1.3)  . 


7 .  Concluding  Remarks 

A  multilevel  approach  for  optimal  design  of  protective 
structures  has  been  presented.  The  solution  method, 
which  is  hased  on  a  simplified  analysis  model  and 
decomposition  of  the  integrated  problem  into  a  number 
of  smaller  subproblems,  is  motivated  by  the  following 
difficulties : 

a.  Nonlinear  dynamic  analysis  is  needed  to  describe 
the  structural  response  even  when  approximate 
models  are  used. 

b.  The  various  types  of  design  variables  are  of 
fundamentally  different  nature  from  both  the 
physical  and  the  mathematical  points  of  view. 

c.  The  objective  function  is  neither  differentiable 
nor  continuous. 

d.  The  problem  size  (numbers  of  variables  and 
constraints)  may  be  large  in  practical  problems. 

A  simple  optimization  model  is  proposed  for  the  cross- 
sectional  variables.  Preopt imi zed  elements  are 
introduced  for  a  set  of  given  element  configurations 
and  loadings.  This  information  is  then  used  for  effi¬ 
cient  optimization  of  the  higher  level  variables.  The 
approach  is  general  and  is  not  restricted  to  a  speci¬ 
fic  problem,  analysis  model,  or  optimization  algorithms. 
Rather,  other  types  of  structure  (such  as  steel  or 
composite  structures  having  non-rectangular  shape), 
different  analysis  models  (dependong  on  the  type  of 
approximations  and  simplifications  used)  and  optimi¬ 
zation  methods  (such  as  optimality  criteria)  can  be 
employed. 

The  numerical  examples  indicate  that  efficient  solu¬ 
tion?  which  do  not  involve  much  computational  effort, 
can  be  achieved  for  complex  optimal  design  problems  by 
the  proposed  approach. 
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